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Using Kubo's linear response theory, we derive expressions for the frequency-dependent electrical 
conductivity (Kubo-Greenwood formula), thermopower, and thermal conductivity in a strongly 
correlated electron system. These are evaluated within ab initio molecular dynamics simulations in 
order to study the thermoelectric transport coefficients in dense liquid hydrogen, especially near the 
nonmetal-to-metal transition region. We also observe significant deviations from the widely used 
Wiedemann-Franz law which is strictly valid only for degenerate systems and give an estimate for 
its valid scope of application towards lower densities. 
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I. INTRODUCTION 

An accurate description of dense hydrogen is of funda- 
mental interest and has wide applications in astrophysics. 
Prominent examples are models of planetary interiors 
which consist predominantly of hydrogen and helium as 
in the case of Jupiter-like planetsP^ Special attention 
has been paid on the high pressure phase diagram of hy- 
drogen and its isotopes. Major problems in this context 
are the slope of the melting line^^ and the transition 
pressure to solid metallic hydrogen,^ which is expected 
near 4 MbaiESl at T = K. Such extreme conditions are 
experimentally still not feasible, yet. 

Interestingly, metallic hydrogen has first been veri- 
fied in the l iquid at about 1.4 Mbar and few thousand 
Kelvin JDJ12I The problem whether or not this transition 
is of first order has been discussed over decades. While 
numerous models within the chemical picture predict al- 
most invariably a pronounced first-order phase transi- 
tion with a critical point at about 10000-15000 K and 
0.5 Mbarp^HIH only few ab initio simulations indicate 
such a behaviorisms] A first-order liquid-liquid phase 
transition with a critical point near 2000 K and 120 GPa 
has been predicted from both quantum Monte Carlo 
(QMC) and finite-temperature density functional theory 
molecular dynamics (FT-DFT-MD) simulations.^ These 
results were confirmed almost simultaneously by exten- 
sive FT-DFT-MD simulations^ which predict a critical 
point at about 1400 K and 132 GPa, i.e. at a somewhat 
lower temperature. A detailed analysis of the changes in 
the structural and electronic properties with density and 
temperature clearly shows that the nonmetal-to-metal 
transition in dense liquid hydrogen drives this first-order 
phase transition. 

The transport coefficients yield valuable information 
on the state of the strongly correlated liquid. For in- 
stance, the conductivity changes drastically along the 
nonmetal-to-metal transition over many orders of magni- 
tude which could already be verified experimentally!^^ 
These measures have the potential to accurately char- 
acterize subtle changes of the electronic structure with 



density and temperature and, especially, alongside the 
liquid-liquid phase transition in hydrogen. We therefore 
calculate the complete set of thermoelectric transport 
coefficients (electrical and thermal conductivity, ther- 
mopower) via FT-DFT-MD simulations for a wide range 
of densities (p = (0.05 — 20) g/cm 3 ) and temperatures 
(T = (1 — 50) x 10 3 K) and pay special attention to the 
nonmetal-to-metal transition region. The focus of our 
work is below 20 g/cm 3 since higher densities are of pri- 
mary importance for the physics of inertial confinement 
fusion and were studied previously!^! 

The transport coefficients as well as their changes with 
density and temperature are again important for appli- 
cations in astrophysics. For instance, a boundary be- 
tween a nonconducting outer and a metallic inner enve- 
lope is usually assumed in interior models of gas giants 
such as Jupiter; its location should in principle be de- 
termined from the nonmetal-to-metal transition in hy- 
drogen but is usually a free parameterPSl Furthermore, 
it has been shown that demixing of helium from hydro- 
gen occurs in H-He mixtures at megabar pressures due to 
the nonm etal-to-metal transition in the hydrogen subsys- 
temPSEIl The treatment of this effect is essential in order 
to explain the excess luminosity of Saturn and its ageP^ 
Knowledge on the transport coefficients of dense liquid 
hydrogen is furthermore required for dynamo simulations 
of planetary magnetic fieldsP^ The dynamo is driven by 
convection of conducting material deep in the interior 
so that the electrical and thermal conductivity are im- 
portant input quantities. Another interesting problem is 
the formation of giant planets out of a protoplanetary 
disc. The conductivities change by many orders of mag- 
nitude during the accretion process as a consequence of 
the density and temperature increase which could affect 
the radiation hydrodynamics of the collapsing disc and 
the luminosity of the young protoplanet strongly. 

There are a number of theoretical models that predict 
the electrical and thermal conductivity with different as- 
sumptions about the el ectron ic and ionic structure and 
their mutual interaction ! 16 ! 30 " 38 ! Especially in a strongly 
coupled system it is difficult to calculate the transport 
properties accurately because the ionization degree and 
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the effective two-particie scattering cross sections are not 
well defined. Strong ion-ion correlations, the dynamic na- 
ture of screening and exchange effects in the electron sys- 
tem as well as quantum effects such as Pauli blocking re- 
quire a consistent quantum-statistical approach. There- 
fore, we calculate the transport coefficients within FT- 
DFT-MD simulations. This method has demonstrated 
its capacity to provide accurate data for the strongly cor- 
related quantum regimePsHUl 

For simple metals the relation between electrical 
and thermal conductivity is described by the famous 
Wiedemann-Franz law^ using a fixed value, the Lorenz 
number Lq, as proportionality constant. The benefit 
from this relation is the ability to obtain the thermal 
conductivity from the electrical conductivity very easily. 
For this procedure the Lorenz number L(T, p) has to be 
known for a given density p and temperature T of the 
system. By calculating the Lorenz number ab initio the 
region where the original Wiedemann-Franz law is valid 
can be identified. 

Our paper is organized as follows. We outline the 
FT-DFT-MD method in Sec. |TT| A generalization of the 
Kubo-Greenwood formula to calculate the Onsager coef- 



ficients Lik is given in Sec. Ill The results for the electri- 
cal and thermal conductivity, the thermopower, and the 
Lorenz number in dense liquid hydrogen are presented in 
Sec. |IV| for a wide range of densities and temperatures. 
In particular, the behavior predicted for the nonmetal- 
to-metal transition region is discussed. Conclusions are 
given in Sec. fVl 



II. THEORETICAL METHOD 



the Baldereschi mean value poi ntP^ w hich proved to yield 
well converged simulation runs |22Mll Several test calcula- 
tions showed that higher efforts are necessary only in the 
vicinity of the phase transition.^ 

The electronic transport coefficients are subsequently 
calculated by evaluating the respective transport formu- 
las, see section below.E31 This is done for 10 to 20 ion 
configurations from the equilibrated MD simulation us- 
ing Monkhorst-Pack k-point meshes^ of 3 x 3 x 3 to 
6 x 6 x 6 to reach the convergence. 

The exchange-correlation functional is the most criti- 
cal input in FT-DFT calculations. Here we use the ap- 
proximation of Perdew, Burke, and ErnzerhoP^ which 
is numerically efficie nt. Thi s functional has been cho- 
sen in similar studies ^ * 57 * 58 ! and reasonable results are 
in general exp ected in the metallic and high-temperature 
plasma regimeP^lSSI In the semiconducting region the ob- 
tained conductivities may be overestimated by using this 
functional but will yet be useful for many practical ap- 
plications. 

Note that the ionic contribution to the transport coef- 
ficients is out of the scope of this work and is therefore 
consistently neglected. 



III. TRANSPORT PROPERTIES AND 
ONSAGER COEFFICIENTS 

In linear response theory (LRT) the response of an 
isotropic single component system of charge carriers 
(electrons with charge q = — e in our case) to an electric 
field E and a temperature gradient VT is expressed by 
the electric current J e and the heat current 3 q through: 



We use the FT-DFT-MD framework which combines 
classical molecular dynamics simulations for the ions with 
a quantum treatment of the electrons based on FT- 
DFT ElHl6] wmcn i s implemented in the VASP 5.2 program 
packageP^l i n the FT-DFT the Coulomb interactions 
between the electrons with the ion s are t reated using the 
projector-augmented wave methocP^U at densities be- 
low 9 g/cm 3 with a converged energy cutoff of 800 eV. 
At higher densities it was necessary to perform all cal- 
culations with the Coulomb potential which required a 
substantial higher cutoff of 3000 eV. 

The FT-DFT algorithm is used to derive the forces 
that act on the ions via the Hcllmann-Feynman theorem 
at each MD step. This procedure is repeatedly performed 
in a cubic simulation box with periodic boundary condi- 
tions for several thousand MD time steps of 0.1 fs to 1 fs 
duration so that the total simulation time amounts up 
to 10 ps. The ion temperature is controlled with a Nose 
thermostat.^! 

Convergence was checked with respect to the particle 
numbers, which vary between 64 and 512 atoms depend- 
ing on the density, the k-point sets used for the eval- 
uation of the Brillouin zone, and the energy cutoff for 
the plane wave basis set. For the simulations we chose 
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The electrical conductivity u, thermal conductivity A, 
and thermopower a are then given by the Onsager coef- 
ficients L mn in the following way: 
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The Lorenz number is the ratio between electrical and 
thermal conductivity divided by the temperature, 



L 



A 



k%Ta 



(4) 



and is, according to the Wiedemann-Franz lawp 21 con- 
stant in the limit of high density, where it can be cal- 
culated as Lo = 7r 2 /3 by means of the Sommerfeld 
expansion.^ 

Within the framework of Kubo's quantum-statistical 
LRT, which is described for instance in Ref. and [Ml 
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the following expressions are obtained for the frequency- 
dependent Onsager coefficients L mn {uS): 



Lmn M = 7^7 (J™ {t - ihr) ; J„)w+j 



(5) 



The current-current correlation functions are given as: 

oo 

(J m (i - ihr) ■ 3 n ) u+lE = lim J dt e^+^ 



p 

J dr Tr |p Jm (t - ihr) ■ J„| 



(6) 



where /3 = (^T)" 1 is the inverse thermal energy. The 
limit e — > has to be taken after the calculation of the 
thermodynamic limit, which is done by evaluating the 
trace. The statistical operator of the equilibrium Qq con- 
tains the Kohn-Sham Hamilton operator iJxs- 

The time-dependent current operators within the 
Heisenberg picture arc defined as: 



J m (t - ikr) = e l(t-^T)^K S j m g-i(t-iftr)H KS 



(7) 



The electric current operator J e = Ji and the heat cur- 
rent operator 3 q — J2 read in second quantization of spin 
degenerate Bloch states: 

J m = J2 (ML|kV)a]> k v' ■ (8) 

Here k is the wave number and v the band index. The 
electric and heat current operators are given as 

Q „ 
— P , 
m e 

1 Hp + pH 
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h e p , 



(9) 
(10) 



where h e is the enthalpy per electron, see Ref. |M] for 
further information on the definition of the currents. 

The eigenvalues of the Hamiltonian can be evaluated 
so that Eq. ([8| can be simplified, 
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(12) 



is used. After inserting Eq. (11) into Eq. ^ and some 
operator algebra the trace can be evaluated according to 
Wick's theorem as: 



Tr {eoaL«k'>PM« P >'} = 

/ki/^k.k 7 fpfi^p,p' ^fJiffi' 

+ /kz^k,p' S v y (1 — fpfj,) <5p,k' $p t i/ f 



(13) 



The first term vanishes and the Fermi functions are de- 
fined as f kl/ = (e^C^-Me) + Altogether the trace 
reads 



4—m—n 



Tr {g J m (t - i hr) ■ J„} = £ e^ R ^ AE 

e kpz^/i 

xfr* (1 - i w ) (HpIpm) • (p/^IpIMC^^L , (14) 

with AE = E kv — E pf2 . Now the r integration can be 
performed: 
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In the same way the second integral can be solved: 
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limy dteH+^+^ AB ) t = 
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(16) 



Only the real part of this equation is considered here, 
since the imaginary part can be calculated more easily 
with a Kramers-Kronig relation. The Onsager coeffi- 
cients now read: 

7i7rg 4 ~ m ~" . e@ AE — 1 
Lmn = 3Vm 2 ^ AE ^ ku _ ^ ^ 

x (k^lplp/x) • (pv\p\kv)6(AE + hu>)eZ- p y p -l. 

In position representation the matrix elements have the 
following form: 

(kv\f»\pfi) = 6 p , k Mik^ + ^ Jd 3 rul v (r)pu kfl (r) 

Here V is the volume of the simulation box and the func- 
tions m£ v (r) and Uk l(i (r) are the Bloch factors. Because 
of the first Kronecker symbol the matrix elements are 
diagonal concerning the wave number, which eliminates 
the p sum. The spin summation leads to an additional 
factor of two. After using a relation between Fermi func- 
tions we finally arrive at the following expression for the 
Onsager coefficients: 



L m n ( w )- 



Wmli 



5~>Hp|k M ) • (MpM (19) 



kz^/j- 



C + C (/k, - A*) S (E kfi - E kv - hu) 



The Onsager coefficients ( 19 1 obey the symmetry re- 
lations L mn {u) = L mn (-Jjaiid L mn (uj) = L nm (w). 
The coefficient L\\ is kno wn as the frequency dependent 
Kubo-Greenwood-formulcP^SI an d has b een widely ap- 
plied in FT-DFT-MD simulations! 39 ' 67 ' 68 ' 
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Similar but not identical frequency dependent formulas 
for L12 and L22 were given by Recoules el alW^ but not 
formally derived in their work. The different formulations 
lead to deviations in the results for the thermopower a 
and in the thermal conductivity A at nonzero frequencies. 
Furthermore, the heat current from Ref . 00] contains the 
chemical potential fj, e instead of the enthalpy h e , see Eq. 
(10 1. However, it can be shown that the term propor- 



tional to the enthalpy per particle h e has no influence on 
the thermal conductivity A in one-component systems. 
As a consequence the numerical results of Refs. [M] and 
HOI could be reproduced by Eq. ( 19 ) in the limit of uj — >• 



and are therefore not questioned by the current work. 
Differences, however, occur in the static results for the 
thermopower a. 

The derivation in this chapter can be easily general- 
ized to the spin dependent form as well as to anisotropic 
systems. 



IV. RESULTS 

A. Electrical conductivity 

Fig. [l] shows isotherms of the electrical conductivity 
over a wide range of densities and temperatures. In pre- 
vious work 41 it could already be shown that these theo- 
retical results agree well with data for the reflectivity and 
conductivity derived from shock-wave experiments. Here 
we concentrate on the general behavior of the theoretical 
curves on a large density and temperature scale and give 
results especially in the vicinity of the nonmetal-to-metal 
transition. 

In general the conductivity rises in the whole range 
studied here with increasing density. On the one hand 
this is the result of pressure ionization which increases the 
amount of conducting electrons and on the other hand, 
even in the case of a fully ionized system, a rising density 
of electrons results in a further increase of the conduc- 
tivity. 

At T — 1000 K the liquid-liquid phase tra nsition which 
was previously reported for dense hydroger l 22 ! 23 ! can be 
identified by the steep increase of the conductivity over 
several orders of magnitude in a small density interval 
between (0.7 — 0.9) g/cm 3 . At temperatures higher than 
1500 K the transition is continuous and the increase 
spreads to a larger range in density. At densities below 
this transition the electrical conductivity increases with 
temperature which is caused by thermal ionization of 
neutral particles. The additional free charges contribute 
to the conductivity. 

At densities above this transition the dependence on 
temperature is inverted: the electrical conductivity de- 
creases with temperature, which is typical for metals. 
Increasing temperature broadens the Fermi function and 
therefore allows additional electron scattering processes 
which reduce their mobility. As a result the conductivity 
is lower with increasing temperature. 



10' 

io e 



L0 J 

1 10 4 
a 

D 

10 2 
10 1 

10° 




10 



50000 K 

30000 K 

20000 K 
15000 K 
10000 K 
■ 5000 K 

2000 K 

1000 K 



0.1 



0.5 



1 



p[g/cm ] 



FIG. 1. (Color online) Electrical conductivity in dense hydro- 
gen for different temperatures: from the region of the liquid- 
liquid phase transition up to the dense plasma. 
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FIG. 2. (Color online) Electrical conductivity for different 
densities versus temperature. 



This behavior is also illustrated in Fig. [2] which shows 
the electrical conductivity along isochores for different 
temperatures. The conductivity decreases for temper- 
atures above 2000 K along the isochores for densities 
higher than 0.9 g/cm 3 which are characteristic of the 
metallic phase. This indicates that most of the system is 
ionized and thus acts metal-like. Looking at lower den- 
sities the conductivity rises along the whole temperature 
range which is due to thermal ionization. The isochores 
clearly show the general behavior of a rising conductivity 
with increasing density. 

Experiments in copper plasmas^ have indicated that 
the electrical conductivity becomes a function of only the 
coupling parameter T for values of T > 10. The plasma 
parameter T is defined by 



(20) 



where n e is the number density of free electrons. Al- 
though such a behavior could not be confirmed later P^l 
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FIG. 3. (Color online) Electrical conductivity as function of 
the coupling parameter Y for different temperatures. The 
plotted data cover a density range of (0.05 — 20) g/cm 3 . 
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FIG. 4. (Color online) Thermal conductivity versus density 
for different temperatures. 



a simple functional form of the electrical conductivity at 
high values of T is still under discussion. To investigate 
whether or not such a simple scaling is valid in dense liq- 
uid hydrogen, the results for the electrical conductivity 
shown in Figs. [T] and [2] are plotted against T in Fig. [3] 
Only for temperatures higher than 10000 K the system is 
strongly ionized. At lower temperatures the occurrence 
of partial ionization prevents a proper calculation of T 
within FT-DFT-MD, because the method does not dis- 
tinguish between bound and free electrons. Therefore we 
do not plot results for lower temperatures in Fig. [3j 

The isotherms appear to be almost parallel in this loga- 
rithmic plot and are clearly separated. Even at the high- 
est available values for r < 60 the isotherms do not tend 
to merge. We conclude that it is not possible to derive a 
simple temperature-independent relation for the conduc- 
tivity that depends solely on the parameter T as it was 
proposed earlier for other metallic liquids. 



B. Thermal conductivity 



This is also shown along the isochores of the thermal 
conductivity that are plotted in Fig. [5j These curves de- 
pict likewise that the thermal conductivity rises invari- 
ably with increasing density and temperature. 
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FIG. 5. (Color online) Thermal conductivity versus temper- 
ature for different densities. 



The isotherms of the thermal conductivity are plotted 
in Fig. [4] They show a similar behavior as the electri- 
cal conductivity and indicate a sharp nonmetal-to-metal 
transition at a temperature below 1500 K. At this tran- 
sition the thermal conductivity increases over several or- 
ders of magnitude in a narrow density range at about 
0.9 g/cm 3 . The increase in the thermal conductivity is 
due to the growing number of delocalized electrons which 
are produced along the phase transition. Above the crit- 
ical temperature this transition becomes broader and is 
caused by a combination of pressure and temperature 
ionization. 

In contrast to the behavior of the electrical conductiv- 
ity, the thermal conductivity does not decrease with tem- 
perature in the metallic phase. The isotherms increase 
systematically with temperature for all densities. 



C. Thermopower 

The thermopower a characterizes the generation of an 
electric field as response to a temperature gradient. For 
most systems this electric field has a direction opposite to 
the temperature gradient which results in a negative ther- 
mopower. The thermopower is most sensitive to changes 
in the electronic structure since it can be expressed as the 
derivative of the logarithm of the electronic conductivity 
with respect to the energy at the Fermi surface.^ Such 
a relation which is also known as Mott formula follows 



from the Kubo-Greenwood equation ( 19 ) in the degener 



ate domain under strong scattering conditions. 

Interestingly, large positive values for the thermopower 
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were measured in fluid mercury^ 72 ! 73 ! near the liquid- vapor 
critical point which is located at T c = 1751 K and 
q c = 5.8 g/cm 3 . In this region, isotherms of the electri- 
cal conductivity near T c show a strong increase with the 
density which is steepest just at the critical density g c . 
This behavior was assumed to be related to fluctuations 
in the electron density which are pronounced near the 
critical point due to critical fluctuations. In particular, 
a zero of the thermopower was observed exactly at the 
critical density. The interesting question arises whether 
or not a zero of the thermopower is a precursor of a first- 
order phase transition in dense liquids which undergo a 
nonmetal-to-metal transition. Wide regions with a posi- 
tive thermopower have been predicted for dense hydrogen 
by a simple chemical modeP^l but were not confirmed in 
an advanced chemical approach^ by ab initio simulations 
or by experiments yet. 
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FIG. 6. (Color online) Thermopower versus density for dif- 
ferent temperatures. 

In Fig. [6] isotherms of the thermopower a are plot- 
ted as function of the density. The symbols repre- 
sent the results from the simulations and the error bars 
show the statistical uncertainties. As guide to the eye, 
polynomial functions were fitted to the numerical re- 
sults. The thermopower is mostly negative and reaches 
a value of about zero at high densities. At lower densi- 
ties the thermopower decreases. The low density limit 
of a — —60.60 /xV/K is known from the Spitzer theory, 
see e.g. Ref. [THl With higher temperatures the neg- 
ative values become systematically larger. We expect 
that these values become smaller again at low densities 
to reach the Spitzer limit. The thermopower shows pos- 
itive mean values below 20000 K and between 0.2 g/cm 3 
and 0.5 g/cm 3 . These values are very small and, in fact, 
about the size of the uncertainties of the respective mean 
values. We therefore cannot conclude definitely that sig- 
nificantly positive values occur within the accuracy of our 
ab initio calculations. To answer the question whether 
or not the zero of the thermopower occurs at the crit- 
ical point of the liquid-liqu id ph ase transition in dense 
hydrogen which is predicte d 22 * 23 ! at about 0.8 g/cm 3 , the 



thermopower has to be evaluated for temperatures below 
2000 K since T c = (1400 - 2000) K. For this region, rea- 
sonable results can be given for the electrical and thermal 
conductivity but not for the thermopower. This is be- 
cause the values of the Onsager coefficients each decrease 
over orders of magnitude in systems with a majority of 
localized electrons while especially the decrease of L i2 is 
two orders of magnitude weaker than that of L\\. This 
behavior on the one hand causes the second contribution 
to the thermal conductivity in Eq. ^ to vanish, because 
Li2 is squared in the numerator. As a result the thermal 
conductivity becomes proportional to L22 which can be 
evaluated successfully in this well as L\\. On 

the other hand this leads to an enormous increase in the 
statistical fluctuations of the thermopower, which causes 
uncertainties that reach hundreds of zxV/K. 
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FIG. 7. (Color online) Thermopower versus temperature for 
different densities. 

In Fig. [7] isochores of the thermopower a are plot- 
ted. The thermopower decreases with higher tempera- 
ture, which is more pronounced the lower the density 
is. Positive values appear at low temperatures, but the 
statistical error is too large to prove this result unam- 
biguously . 



D. Lorenz number 

The Lorenz number Eq. Q describes mainly the re- 
lation between thermal and electrical conductivity. For 
simple metals this relation is described via a constant and 
is known as the Wiedemann-Franz law. For high degen- 
eracy this constant is Lq = 7r 2 /3pU This relation can be 
used to easily obtain the thermal conductivity for metals 
if the electrical conductivity is known. Here we calculate 
the Lorenz number for dense hydrogen in order to iden- 
tify the region where the Wiedemann-Franz law is valid. 
High degeneracy occurs only at sufficiently high densities 
so that we expect deviations from the Wiedemann-Franz 
law at lower densities and higher temperatures. 

The Lorenz number is shown in Fig. [8] for several tern- 
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FIG. 8. (Color online) Lorenz number versus density for dif- 
ferent temperatures. The dashed black line displays the lim- 
iting value of Lo = 7r 2 /3. 

peratures. At densities above 1 g/cm 3 the Lorenz num- 
ber is almost constant and shows no temperature depen- 
dence which indicates that the Wiedemann-Franz law is 
valid there. The value of Lo = 7r 2 /3, which is shown by 
a dashed black line, could be reproduced for each tem- 
perature at densities higher than (1-2) g/cm 3 within the 
statistical uncertainties. This behavior is consistent with 
the metallic-like properties observed already for the con- 
ductivity at high densities. 

At densities below (1-2) g/cm 3 the Lorenz number 
shows strong deviations from the Wiedemann-Franz law 
with a pronounced temperature dependence. In partic- 
ular, the Lorenz number rises strongly with decreasing 
density. This behavior becomes more pronounced at 
lower temperatures. On the other hand, at the lower 
temperatures the validity of the Wiedemann-Franz law 
extends to smaller densities, e.g. down to 0.2 g/cm 3 at 
10000 K and 15000 K. 

Along the isotherms of 30000 K and 50000 K a max- 
imum can be identified. Thus, at lower densities the 
Lorenz number decreases again. At 30000 K this maxi- 
mum is found at a smaller density than at 50000 K and 
shows a larger value. This indicates that with decreas- 
ing temperature the maximum shifts systematically to 
smaller densities and increases. For temperatures below 
30000 K the maximum is not displayed. It appears to 
be situated at densities below 0.05 g/cm 3 , which is the 
minimal density calculated here and which displays the 
current limit of our computational capabilities. However, 
we expect a maximum to appear for all temperatures be- 



cause the Lorenz number is also known exactly at low 
densities from the Spitzer theory for fully ionized plas- 
masPSl The latter predicts a constant low-density value 
of 1.5966, see e.g. Ref.QIl 

Our results indicate that the deviations from the 
Wiedemann-Franz law might reach a full order of mag- 
nitude in certain regions of density and temperature. 
Further research will be necessary to investigate this im- 
portant aspect in more detail. We also expect that cer- 
tain deviations from the Wiedemann-Franz law occur in 
warm dense matter of arbitrary composition. In hydro- 
gen, we predict that the Wiedemann-Franz is valid at 
densities above 2 g/cm 3 for temperatures below 50000 K, 
see Fig. [8] 

V. CONCLUSIONS 

We derived formulas for the Onsager transport coeffi- 
cients L mn within the Kubo theory for application in ab 
initio simulations. Using these expressions we calculated 
the electrical and thermal conductivity as well as the 
thermopower of liquid hydrogen for a wide range of tem- 
peratures and densities in the megabar range. In partic- 
ular, we characterize the nonmetal-to-metal transition in 
hydrogen by observing a rapid increase in both the elec- 
trical and the thermal conductivity. The thermopower 
shows a trend towards positive values in a region where 
the critical point of the liquid-liquid phase transition is 
expected, similar to the behavior of liquid mercury. At 
low temperatures more accurate calculations will be nec- 
essary for this transport coefficient which are beyond the 
scope of the currently available computer capacity. The 
ab initio calculation of the Lorenz number shows in ad- 
dition that the validity of the original Wiedemann-Franz 
law is limited to the metallic-like regime of hydrogen. 
The theoretical framework given here can be applied to 
bulk material calculations for arbitrary materials within 
FT-DFT-MD. 
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